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Realistic 3D radiative MHD simulations open new perspectives for understanding the turbulent dynamics of the solar surface, its coupling to the atmosphere, and the physical mechanisms of generation and transport of non-thermal energy. Traditionally, plasma eruptions and wave phenomena in the solar atmosphere are modeled by 
prescribing artificial driving mechanisms using magnetic or gas pressure forces that might arise from magnetic field emergence or reconnection instabilities. In contrast, our 'ab initio' simulations provide a realistic description of solar dynamics naturally driven by solar energy flow. By simulating the upper convection zone and the solar 
atmosphere, we can investigate in detail the physical processes of turbulent magnetoconvection, generation and amplification of magnetic fields, excitation of MHD waves, and plasma eruptions. We present recent simulation results of the multi-scale dynamics of quiet-Sun regions, and energetic effects in the atmosphere and compare 
them with observations. For the comparisons we calculate synthetic spectro-polarimetric data to model observational data from SDO, Hinode, and the BBSO New Solar Telescope. 
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3D rectangular geometry 

Fully conservative, Fully compressible 

Fully coupled radiation solver: 

LTE using 4 opacity-distribution-function bins 
Ray-tracing transport by Feautrier method 
1 8 ray (2 vertical, 1 6 slanted) 
angular quadrature 
Non-ideal (tabular) EOS 
4th order Pade spatial discretization 
4th order Runge-Kutta time integration 
LES-Eddy Simulation options (turbulence models): 
Compressible Smagorinsky model 

Compressible Dynamics Smagorinsky mode 
(Germano et a., 1 991 ; Moin et al, 1 991 ) 

MHD subgrid models (Balarac et al., 2010) 
DNS+Hyperviscosity approach 
MPI parallelization (plane and pencil versions) 


Magnetized vortex tubes and flow eruptions 

a) 1 b) 


B z0 =10G, X 0 =6301.5A, -0.06A, At=10s 
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Panel (a) shows a vertical cut through a vortex tube which drives an eruption and generates shock waves in the solar 
chromosphere. Horizontal cuts (panel (b) show a fraction of the domain at the photosphere (top row) and at 625 km above the 
solar surface (bottom row) for different quantities (left to right): vertical velocity, density, temperature, and vertical magnetic field. 
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Basic equations 

The equations are the grid-cell averaged 
Conservation of mass: dp 
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Conservation 
of momentum: 
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Comparison of the contributions to 
the flow ejection of the hydrodynamic 
(red curves) and magnetic (Lorentz 
force, blue curves) vertical 
accelerations (normalized by the 
gravity acceleration) at different 
altitudes: (a) h = 780 km, (b) h = 

300 km, (c) h = 0 km (photosphere), 
and (d) h = -240 km. The dashed 
curves correspond to the vertical c) 12 ^ 
velocity, V z , given for reference at the 
same layers. 
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z = 300 km i 


Small-scale eruption signatures reflected in Stokes I, V and (U+Q). Panels show zoomed images of l/lc, V/lc and (U+Q)/lc for A?^=-0.06A (X 0 =6301 .5A). 
The time-difference between snapshots is 10s. White circle shows the location of a probe for studying the evolution of the Stokes profile shape. 

Modeling center-to-limb variations of the HMI observables 
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Conservation of magnetic flux 


d B t / 
d t 


Structure of 
B j B j the flow 
during 
eruption. 


+ 


u B - uB - 


c 


j i 


i j 


\ 


4 JTO 


(«. - « ) 


‘SolarBox’ code 

1) 80x80x26 Mm 

2) 6.4 x 6.4 x 6 Mm 


Initial 

magnetic 

field: 
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Local dynamo simulations 
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Horizontal snapshot of a fraction of the computational domain in the 
photosphere for: a) the vertical velocity; b) magnetic field strength; c) 
vertical component of the electric current density; d) logarithm of the 
turbulent magnetic Prandtl number; e) turbulent magnetic diffusivity; 
and f) time-derivative of the magnetic energy density. 



Temporal profiles of 
variations of the mean 
vertical velocity at 
different levels below o.e 
the surface and in the 
atmosphere. The thick 
black curve shows the 
variations in the 
photosphere layer. 
The height difference 
between the curves is 
60 km. 
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Square-root of the power spectral density of the 
simulated Doppler-shift signal for the viewing 
angles. The spectra reveal the line asymmetry, 
which is important for helioseismology and ring- 
diagram analyses. 

Panel (a): Doppler-shift velocity variations as a 
function of time. Panel (b): time-averaged Doppler - 40 o 

shift velocity as a function of the viewing angle. t, min lag time, sec 

Panel (c): Doppler-shift velocity fluctuations after removing the angular dependence. Panel (d): cross-correlations of the 
oscillations observed at different angles relative to the disk center. Different line colors correspond to different angular 
distances from the disk center: 0° (blue curves), 30° (green), 45° (orange) and 60° (red curves). 



Topology of the magnetic field lines above the photosphere in 
the local dynamo simulations. The horizontal plane shows the 
distribution of the vertical magnetic field in the photosphere. 
Red color corresponds to positive polarity, blue color to 
negative polarity of the vertical magnetic field. The range of 
field strength is from 800 G to 300 G. 



Mean vertical profiles of the unsigned magnetic field 
components (averaged over 2 hours): vertical component 
(red curve) and transverse component (blue). Error bars 
correspond to 0.5a. 
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Ratio of the mean vertical and transverse components 
magnetic fields as function of depth below the photosphere. 
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Comparison of Stokes l/lc and V/lc for four Flow eruption reflected in Stokes I and V. Different 
lines during the flow eruption. Time and panels illustrate the evolution of Stokes profiles l/lc and 
location of the region of sample is indicated V/lc for X 0 =5250, 6301, 6173 anA. Red-line profiles 

correspond to time and location marked by circle above. 
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Different stages of a flow ejection. 
Black streamlines illustrate the 
velocity field in the vicinity of a vortex 
tube. Semitransparent light gray 
surface corresponds to a constant 
temperature of 6400 K. Yellow and 
blue isosurfaces show the vertical 
pressure gradient divided by density 
(-Vp/p) is equal to 5-1 0 4 cm/s 2 
(yellow color) and -5-1 0 4 cm/s 2 (blue 
color). 
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Comparison of the 
synthetic line profiles 
obtained from our 
simulations for two iron 
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Ridges in nu, / - plane calculated from the synthetic data 

(Fe I, 6173A): Doppler shift (top) and continuum intensity e 
( bottom) for 0° and 60° from the solar disk center. 
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abundances: 7.43 and l 20 ^ 


7.50 dex, and from 
observations with 
Fourier Transform 
Spectrometer (FTS 
solar atlas) for the line 
61 73.3 A. 
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Variation of the height of 
formation of the 6173.3 A 
line as a function of the 
angular distance from the 
: solar disc center for the 
full-resolution (red dots) 
and HMI resolution 
(blue). Both sets of data 
1 are fit with a second 
order polynomial. 
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Doppler-shift snapshots calculated from the simulated profiles of the 
Fe I, 61 73A line for the solar disc center for models with full numerical 
resolution (12.5 km), and with HMI resolution (-315 km/pixel). 


Conclusions 

Realistic 3D radiative hydrodynamic simulations allow us to investigate important links between solar 
observations and the physical and dynamical properties of solar conditions. This is achieved by taking into 
account spectral line formation conditions, line-of-sight effects, and instrumental spatial and spectral resolutions. 
It is important to model the observational and data analysis procedures and calculate 'observables' using the 
simulation data to achieve direct comparison of observations and numerical models . 

The small-scale dynamics are not easily observed in broad-band images even with high observational spatial 
resolution. We examine four lines (5250A, 6173A, 6301 A and 15648A) for studying the temporal evolution of the 
Stokes profiles during small-scale flow eruptions in magnetic field regions. The synthetic Stokes profiles reveal 
interesting variations in the form of ‘running-wave' perturbations across the line with time - reflecting the 
bidirectional flow dynamics of the eruptions. The Stokes images show ubiquitous small-scale elongated structures 
which originate in the photosphere. Some of these can be identified as ultrafine loops above the granulation. 

We found that the primary topology of the dynamo-generated magnetic field is represented by compact magnetic 
loops appearing as bipolar structures in the intergranular lanes, where the height of the loops is greater for 
stronger magnetic field concentrations at the foot-points. The formation of magnetic loops in the solar atmosphere 
causes the observed anisotropy as a function of height of the vertical and transverse small-scale magnetic fields. 

















